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Abstract 


Proper water management in direct methanol fuel cells (DMFCs) is very critical and complicated because of many interacting physicochemical 
phenomena. Among these, the liquid saturation in the cathode side is believed to have a very strong effect on water crossover through the membrane, 
a key parameter to determine water balance between the anode and cathode. In this paper, based on an interfacial liquid coverage model implemented 
in a three-dimensional (3D) two-phase DMFC model, the liquid saturation variations in the cathode are examined in detail and their effects on the 


net water transport coefficient through the membrane discussed. 
© 2006 Elsevier B.V. All rights reserved. 
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1. Introduction 


The direct methanol fuel cell (DMFC) is considered as a 
leading contender for next-generation portable and micro power 
sources, offering a combination of simplicity, robustness and 
high energy density. Recent research indicates that water man- 
agement [1,2] is one of most critical challenges for DMFCs to 
compete with lithium-ion batteries. 

Constrained by the methanol crossover problem, the anode 
fuel solution has been very dilute, meaning that a large amount of 
water needs to be carried in the system and therefore reducing the 
energy content of the fuel mixture. In addition, for each methanol 
molecular reacted in the anode catalyst, 16 water molecules 
will be lost from the anode, assuming that the electro-osmotic 
drag coefficient of water is equal to 2.5 per proton through a 
thick membrane, such as Nafion 117 [3]. On the other hand, 
inside the cathode, there are 15 water molecules transported 
from the anode and 3 additional water molecules produced by 
consuming 6 protons generated from oxidation of 1 methanol. 
How to minimize water removal from the cathode and subse- 
quent recovery externally to replenish the anode without caus- 
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ing severe cathode flooding becomes an important engineering 
issue. 

Formally, the net water flux through the membrane, caused 
by diffusion, electro-osmosis, and hydraulic pressure differential 
across the membrane, can be quantified in terms of the net water 
transport coefficient, expressed as a = NIRO / NE [4,5]. The 
ideal water transport condition, i.e. a water neutral cell, can be 
achieved when «œ = — 1/6 [2]. In other words, no water is needed 
from the anode feed and one water molecule needed to oxidize 
methanol comes from the product water generated on the cath- 
ode side. The magnitude and spatial distribution of œ are thus 
very important for the design of innovative water management 
strategies in DMFCs. 

Based on the theory of liquid water transport in proton 
exchange membrane fuel cells (PEMFCs) [6,7], a unique MEA 
structure has been developed, which utilizes a thin membrane 
(e.g. Nafion 112) to promote water back flow from the cathode 
to anode. Such MEAs exhibit extraordinarily low œ and hence 
are generally termed low-a MEA technology [1,2,5]. 

Recent visualization data from H2/air PEMFCs show that 
in the high current density and/or low gas flowrate conditions, 
water produced from the oxygen reduction reaction is removed 
from the GDL and gas channel in the liquid phase. Liquid water 
transport through hydrophobic GDL is driven by the gradient in 
capillary forces. At the hydrophobic surface of the GDL facing 
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Fig. 1. Image of the cathode backing layer covered by water droplets in seven- 
channel flow field with a total active area of 3 cm? [8]. 


the gas channel, liquid water emerges in the form of droplets, 
as shown in the experimental image of Zhang et al. [8], which 
is reproduced here in Fig. 1. The droplet detachment diameter 
along with the droplet population density at the GDL surface will 
result in an interfacial liquid saturation, which is a key parameter 
to determine the level of flooding inside the GDL and catalyst 
layer. This same concept of interfacial coverage is believed to be 
important in a portable DMFC because much water exists in the 
DMFC cathode due to water crossover through the membrane. 

In tandem with experimental water transport research, var- 
ious mathematical models were also developed in order to 
achieve a deeper understanding of water behavior in the PEFC 
[6,7,9]. However, mathematical modeling of water transport in 
DMFCs has received little attention if any, due to the recog- 
nization of its importance only most recently. For the first time, 
Liu and Wang [10] presented a 3D, two-phase DMFC model, 
which includes flow channels, backing and catalyst layers on 
both anode and cathode sides, and the membrane as a single 
simulation domain, that elaborates water transport. This model 
is extended to further study the net water transport coefficient 
distribution and interfacial liquid water coverage effect in the 
present paper. A more complete review of DMFC modeling in 
the literature is recently given by Wang [11]. 

In the following, a 3D, two-phase mathematical model of 
DMEFCs with the above-mentioned water transport features will 
be presented, and then solved by a commercial flow solver, 
Fluent®. Model results will be presented to illustrate the net 
water transport coefficient distribution in a 3D geometry and 
water transport behavior under different liquid water coverage 
conditions. 


2. Physical model 


The two-phase steady-state DMFC model in this work con- 
sists of non-linear, coupled partial differential equations of con- 
servation of mass, momentum, species and charge. These can 


be presented in concise form as follows [10]: 
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where p, u, c and e denote the two-phase mixture density, 
velocity vector, molar concentration of species k and electrolyte 
phase potential, respectively. The species considered here are 
methanol, water and oxygen. The various two-phase proper- 
ties, source terms, electro-chemical kinetics and thermophysical 
properties identified for various regions of a DMFC, as well 
as necessary boundary conditions, have been detailed in refs. 
[10,12—14] and thus are not repeated here. 

Liquid saturation is a key parameter in the two-phase flow 
model. Here, we can obtain the liquid saturation from the mixture 
water molar concentration via 


_ cO 
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Zhang et al. [8,15] most recently found that the interfacial liq- 
uid saturation correlates strongly with the channel gas velocity as 
the droplet detachment diameter from the backing layer surface 
decreases with the gas velocity as a result of the balance between 
the drag force of air core flow on the drop and the drop’s sur- 
face adhesion with the cathode backing layer. In addition, since 
liquid tends to form a droplet or a film at the hydrophobic and 
hydrophilic porous surfaces, respectively, the interfacial liquid 
saturation at the cathode backing layer and channel interface 
is also a function of the cathode backing layer wettability, i.e. 
contact angle. Furthermore, the droplet population density at the 
cathode backing layer surface, i.e. the number of sites activated 
with droplet emergence, is related to the local current density. 
Therefore, the interfacial liquid saturation is also a function of 
the current density. To summarize, we have 


Sint = s(u, be, i) (6) 


The exact form of this interfacial saturation function is to be 
determined by experimental measurements for a specific cathode 
porous material system. In the present 3D model focusing on 
exploring fundamental physics and numerical behaviors of two- 
phase dynamics with interfacial liquid coverage at the cathode 
backing layer surface, we implemented a constant interfacial 
liquid saturation for the purpose of a parametric study in this 
work. 
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Table 1 

Three-dimensional cell geometry and operating conditions 

Cell length (m) 0.1 

Cell width (m) 2x 10-3 

Anode channel width (m) 1x 103 

Anode backing thickness (m) 300 x 1076 

Anode catalyst thickness (m) 10 x 1076 

Membrane thickness (m) 50 x 1076 

Cathode catalyst thickness (m) 10 x 1076 

Cathode backing thickness (m) 300 x 1076 x 7 = À 

Cathode channel width (m) 1x 1073 te N 

Operating temperature (°C) 60 

Code ohanielpreadure a) ) 

Flow stoichiometry of anode channel 2 

Flow stoichiometry of cathode channel 3 

Inlet methanol concentration at anode 2000 (2 M) 
(mol m~3) 

Inlet oxygen concentration at cathode 7.68 Fig. 2. Illustration of 3D simulation domain and several planes for presentation 
(mol m~?) of simulation results. 

Inlet liquid saturation at anode (%) 100 

Inlet liquid saturation at cathode (%) 0 (fully humidified air) ation in the cathode catalyst layer and water transport from the 

eas ey the 0, 5, 10,20 anode to cathode, the highest liquid saturation is about 7% in 

Cell operating voltage (V) 0.4 the cathode backing layer under given operating conditions. In 


3. Results and discussion 


As discussed, water transport through the membrane is gov- 
erned by electro-osmotic drag, diffusion and hydraulic perme- 
ation. The electro-osmotic drag coefficient is only a function 
of temperature and thus a constant at any fixed operating cell 
temperature at steady state. Water diffusion and hydraulic per- 
meation through the membrane are highly dependant on the 
liquid saturation of both anode and cathode sides. With the liquid 
coverage model applied at the interface of the cathode backing 
layer and flow channel, water diffusion driven by the water con- 
tent gradient and water permeation by the hydraulic pressure 
gradient vary significantly and therefore, affect the net water 
transport coefficient, œ, considerably. 

The predictions of water crossover coefficient using a 1D ver- 
sion of the present model has been found to be in good agreement 
with experimental data [16], demonstrating the validity of this 
physical model. In the following, focus is placed on discussing 
3D water transport simulation results for a 3D geometry and 
operating conditions described in Table 1. To clearly understand 
numerical results shown in this section, several representative 
planes in the 3D domain are described in Fig. 2 to facilitate the 
presentation of numerical results. 

The liquid saturation distribution in the cathode backing 
layer without interfacial coverage is shown in several planes 
in Figs. 3—5. In this case, the liquid saturation distribution in 
the backing layer is greatly affected by the land of the current 
collector. The liquid saturation in the channel area of the back- 
ing layer is lower than that under the land. Since liquid water in 
the backing layer must be transported from the land area to the 
channel area for removal by channel gas, the liquid saturation 
in the land area must be higher in order to provide the capillary 
forces to this lateral liquid water motion [14]. With water gener- 


Fig. 4, it is evident that the gradient of liquid saturation in the 
through-plane direction is much stronger than that in the flow 
direction. And, in the Plane X3, parallel to the membrane, the 
liquid saturation is lower in the channel region but rising around 
the edges due to the effect of the current collector lands. 

The net water transport coefficient contour without liquid 
coverage is shown in Fig. 6. In contrast with the 1D model 
assumption, the net water transport coefficient is found to be 
quite non-uniform, varying from 1.2 to —1.0. This indicates 
that net water transport through the membrane in the two- 
dimensional plane is not only from the anode to cathode, but 
also from the cathode to anode, depending on the location, local 
liquid saturation and current density. In the large middle region 
of the membrane, the net water transport coefficient is relatively 
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Fig. 3. Liquid saturation in the cathode backing layer of Plane Y, without liquid 
coverage. 
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Fig. 4. Liquid saturation in the cathode backing layer of Plane Z, without liquid 
coverage. 


uniform and almost constant, while it shows a large gradient and 
becomes even negative in the regions covered by lands. 

To further investigate the reason for such a water trans- 
port coefficient distribution, water transport through the mem- 
brane due to diffusion and hydraulic permeation is shown in 
Fig. 7. Here, the water transport coefficients due to diffusion 
and hydraulic permeation are defined as œqiff = Ne! NH 
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Fig. 5. Liquid saturation in cathode backing layer (Plane X3), without liquid 
coverage. 


and op) = N /N HY respectively. Therefore, the net water 


transport coefficient can be written as [10] 


H20 


a= ng + Qdiff — Op! (7) 


where no is the water electro-osmotic drag coefficient through 


the membrane, which is mainly depending on the operating 
temperature. It is clear that both water diffusion and hydraulic 
permeation have very uniform distribution in the most middle 
region of the membrane controlled by rather uniform current 
density distribution. Due to the stronger effect of hydraulic per- 
meation from the cathode to anode through the membrane in 
the areas covered by lands, net water transport is pointing from 
cathode to anode, that is, the net water transport coefficient is 
negative. 

Under the same operating conditions, the liquid saturation 
distribution in the cathode backing layer with 10% liquid cover- 
age is shown in Figs. 8 and 9. Compared with no coverage results, 
the liquid saturation in the cathode backing layer is almost con- 
stant, pre-dominated by the interfacial liquid coverage of 10%. 
The almost constant distribution is due to the relatively large 
capillary diffusivity with about 10% liquid saturations in the 
backing layer [11]. With such a large capillary diffusivity, the 
liquid saturation distribution is fairly uniform in Plane X3 in the 
cathode backing layer, shown in Fig. 9. In other words, the land 
effect becomes negligible for the cathode liquid saturation with 
the interfacial liquid coverage above 10%. 

The net water transport coefficient distribution with 10% liq- 
uid coverage is shown in Fig. 10. The general trend is the same as 
that without liquid coverage, however, with lower coefficients 
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Fig. 6. Net water transport coefficient, a, distribution through the membrane, 
without liquid coverage. 
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(a) By diffusion (b) By hydraulic permeation 
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Fig. 7. Water crossover due to diffusion (a) and hydraulic permeation (b) through the membrane, without liquid coverage. 


in the entire membrane region. The lower net water transport 
coefficient can be explained by stronger hydraulic permeation 
in Fig. 11(b). For both cases, water diffusion distribution from 
the anode to cathode is about same, while average hydraulic per- 
meation through the membrane with 10% coverage is about 20% 
higher than no coverage case. Therefore, the net water transport 
coefficient becomes smaller with 10% liquid coverage than no 
coverage. 
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Fig. 8. Liquid saturation in the cathode backing layer of Plane Z, with liquid Fig. 9. Liquid saturation in the cathode backing layer (Plane X3), with liquid 


coverage of 10%. coverage of 10%. 
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Fig. 10. The net water transport coefficient, œ, distribution through the mem- 
brane, with liquid coverage of 10%. 


1.15 
1.11 
1.07 
1.03 
0.99 
0.95 
0.91 
0.87 
0.83 
0.79 
0.76 
0.72 
0.68 
0.64 
0.60 
0.56 
0.52 
0.48 
0.44 
0.40 


(a) By diffusion 


1 70.155 
E 0.9 E Net water transport coefficient E 
S 0.8 =*= Average current density +0.15 z 
E = 
ge | 0.145 È 
t 06 i= 
Q Š 
2 0.5 +0.14 © 
S 
© 04 g 
ee | 5 
5 03 0.135 3 
8 © 
igs +013 $ 
© 01 ` o 
a é 

0 + + + 0.125 


$ 
0 5% 10% 20% 
Interfacial liquid coverage 


Fig. 12. Net water transport coefficients and average current densities with var- 
ious liquid coverage saturations. 


Using different liquid coverage saturations, the net water 
transport coefficients and average current densities are predicted 
in Fig. 12 by the 3D model. It is evident that the net water 
transport coefficient is highly dependant on the liquid coverage 
condition at the cathode backing layer surface, varying from 
0.94 under the no coverage condition to 0.22 under the liquid 
coverage of 20%. The average current density shows relatively 
small variation with higher liquid coverage saturations under 
the given operating conditions. However, it may possibly show 
significant effect when oxygen transport is critical, for exam- 
ple, under high current density conditions or with air-breathing 
cathodes in portable DMFCs. 
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Fig. 11. Water crossover due to diffusion (a) and hydraulic permeation (b) through the membrane with liquid coverage of 10%. 


W. Liu, C.-Y. Wang / Journal of Power Sources 164 (2007) 189-195 195 


4. Conclusions 


In this paper, an interfacial liquid coverage model applied at 
the interface between the cathode backing layer and flow chan- 
nel is numerically investigated based on a 3D two-phase model, 
and its effects on the net water transport coefficient distribu- 
tion in a DMFC are explored under typical operating conditions 
for portable applications. It can be seen that interfacial liquid 
coverage has a profound effect on the net water transport coef- 
ficient through the membrane by affecting water diffusion and 
hydraulic permeation. With an accurate correlation of liquid cov- 
erage from experimental measurements, water balance between 
the anode and cathode in the DMFC can be effectively tailored to 
be suitable for use of high concentration fuel, without sacrificing 
cell performance. 
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